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ABSTRACT 


A  unified  approach  is  applied  to  the  derivation  of  a  number  of  formulas 
for  the  probability  of  signal  detection  and  the  probability  of  false  alarm. 

The  context  is  incoherent  integration,  with  fluctuating  signal-to-noise  ratios 
and/or  fluctuating  thresholds.  The  standard  results  are  obtained  and  extended 
to  more  general  fluctuation  models.  A  fundamental  duality  is  established 
between  fluctuating  signals  and  fluctuating  thresholds  and  used  to  simplify 
the  derivations.  Also  included  is  an  expression  of  the  cumulative 
F-distribution  as  a  finite  sum  of  llarcum  Q-functions. 
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I.  INTRODUCTION 


In  a  wide  variety  of  signal  detection  and  binary  decision  problems,  the 
evaluation  of  performance  reduces  to  the  computation  of  a  number,  P,  defined 
as  the  probability  that  one  random  variable,  u.,  exceeds  another,  called  y. 

In  symbols, 

P  =  Frob  {u  >  y}  , 


and  u  represents  a  processor  output  which  is  being  tested  for  signal  preserve, 
while  y  is  a  threshold,  which  may  be  a  constant.  In  all  cases  considered 
here,  u  and  y  are  independant,  and  u  is  the  result  of  incoherent  summation  of 
N  complex -samples  each  containing  signal  and  noise  components.  These 
quantities,  z^,  represent  I  and  Q  output  samples  of  some  coherent  processor 
and  they  are  modelled  as  sums  of  signal  and  noise  terms: 

z  =  s  +  w  . 
n  n  n 

The  noise  components  are  Gaussian,  independent  and  with  zero-mean,  all  sharing 
the  same  variance: 

E|w  I2  =  2o2  , 

n 

and  all  having  the  circular  property 

Ew  2  =  0. 
n 


Various  models  will  be  used  for  the  signal  components,  corresponding  to  (and 
extending)  most  of  the  "fluctuating  target"  models  of  the  radar  literature. 
The  random  variable  u  is  the  normalized  sum 


l  N 

u  =  — s-  f  Iz  | 
2„2  nil  n 


(1-1) 
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hence  its  probability  distribution  function  (pdf),  conditioned  on  the  signal 
sequence,  {sn)  ,  is 


fN(u,a)  =  e  U  a(u/a)  1^ 


(1-2) 


In  this  formula. 


a 


N 


l 

n~l 


(1-3) 


and  In  is  the  Bessel  function  of  imaginary  argument.  Since  the  signal 
components  enter  this  conditioned  pdf  only  through  the  sum  a,  the  "random 
signal"  models  will  be  represented  by  postulated  pdf’s  for  a. 

In  a  similar  way,  "fluctuating  thresholds,"  such  as  those  arising  in  CFAR 
problems,  are  represented  by  putting  y=cx,  where  c  is  a  scalar  multiplier, 
and  x  is  a  random  variable.  Various  pdf's  are  postulated  for  x,  including  the 
pdf  f^(x,d),  with  fixed  or  fluctuating  parameter,  d.  Ordinary  "linear  CFAR” 
is  represented  by  this  model  with  d=0. 

'2  3) 

When  a  and  y  are  constant,  P  is  the  standard  Marcum  Q-f unction^  ’  '  which 
must  be  computed  from  an  infinite  series,  used  together  with  a  bound  on 

(4) 

truncation  error.  An  efficient  algorithm  has  been  described  by  Shnidman, 
which  involves  recursive  computation  of  the  terms,  along  with  several 
refinements  which  improve  computational  efficiency.  The  standard  series, 
(Fehlner's  formula,  see  Section  6)  was  used  by  Mitchell  and  Walker as  the 
basis  of  a  systematic  derivation  of  expressions  for  some  of  the  cases 
presented  here.  Because  they  start  with  an  infinite  series  for  the  Q-function 
and  obtain  their  results  by  averaging  over  various  random  models  for  a  and/or 
y,  all  the  expressions  obtained  are  in  infinite  series  form.  Although  all 
these  series  solutions  lend  themselves  readily  to  recursive  computation,  it 
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happens  chat  in  most  of  the  cases  (except  the  original  one,  where  both  a  and  y 
are  constant),  finite-sum  expressions  exist  and  many  of  these  can  be  found  in 
the  literature. ^ 

In  this  study,  we  start  with  an  integral  representation  for  the 
Q-function,  instead  of  an  infinite  series.  Then,  in  analogy  to  the  method  of 
Mitchell  and  Walker,  we  derive  the  desired  results  for  random  a  and/or  y  by 
averaging  this  expression  under  the  sign  of  integration.  Then,  by  simple 
changes  of  the  variable  of  integration,  the  integral  representation  is  itself 
used  to  obtain  the  desired  finite-sum  expressions.  In  all  but  two  cases,  the 
terms  of  these  sums  are  elementary  functions,  easily  computed  recursively.  In 
the  two  exceptions,  the  terms  of  the  sums  are  themselves  Q-functions.  One  of 
these,  actually  representing  the  cumulative  F-distribution,  is  the  source  of  a 
series  of  generalizations  of  the  standard  fluctuation  models. 

When  a=0  in  formula  (1-2),  the  resulting  pdf  will  be  called  g^(u): 


M  JL 

gN(u)  «  fN(u,0)  e 


0-4) 


This  function,  the  chi-squared  with  2N  degrees  of  freedom,  is  used  to 
model  signal  and  threshold  fluctuations  for  the  standard  formulas  derived  in 
Section  4.  Only  integral  values  of  N  are  used,  since  non-integral  values  (as 
in  the  Weinstock  cases^^)  do  not  lead  to  detection  formulas  of  finite  form. 
Fluctuating  signals  are  modelled  by  postulating  that  the  pdf  of  a  is 


Va>  *  I  %<t> 


(1-5) 


where  M  and  b  are  fluctuation  parameters.  In  this  model,  the  mean  value  of  a 


a  =  Mb, 

and  according  to  (1-3), 


(1-6) 


1  N  7 

I  =  -V  l  <s  i 

2c2  n-1  " 


(1-7) 


We  denote  the  average  signal-to*-noise  ratio  for  a  single  sample  by  “SNR,”  so 
that 


SNR  =  | 


for  a  fixed  signal,  and 

SNR  =  |  =  |  b  (1-8) 

for  a  signal  fluctuating  according  to  (1-5).  The  standard  “Swerling  cases” 
correspond  to  M-values,  as  follows:'2^ 


Swerling  Case  M  Value 

1  1 

2  N 

3  2 

4  2N 


When  (1-4)  is  used  to  describe  a  fluctuating  threshold,  we  write 
y  =  cx, 

and  use 

gL(x) 

for  the  pdf  of  x.  The  pdf  of  y  is  then 


f  (y)  =  — 

o  c 


*L<t> 


(1-9) 


4 


This  model  represents  linear  CFAR,  and  the  variable  x  is  the  result  of 

incoherent  summation  of  L  complex  samples  containing  only  noise  of  the  same 

variance  as  that  of  “he  signal  samples.  The  terms  of  this  sum  are  normalized, 

2 

just  as  in  (1-1),  by  the  true  noise  variance,  a  ,  but  these  normalizing 
factors  cancel  when  the  probability  that  u  exceeds  y  is  computed.  If  the 
noise  samples  are  known  to  have  a  different  variance  than  the  signal  samples, 
the  variance  of  the  signal  samples  is  used  for  normalization  and  the  variance 
ratio  is  absorbed  into  the  constant,  c.  Note  that  x  is  the  sum  of  the  noise 
samples,  not  the  average,  which  affects  only  the  significance  of  c. 

The  fluctuation  models  are  generalized  in  Section  5.  First,  the  pdf, 

£o(a>  -H(t- d)-  w-io) 

is  used  to  describe  signal  fluctuations.  The  expected  value  of  a  in  this  case 
is  a  =  (M4d)b,  hence 

SNR  =  |  b(l  +  £)  .  (1-11) 

This  case  is  of  rather  academic  interest,  and  the  dual  problem  with  fixed 
signal  and  threshold  y=cx,  is  more  interesting.  The  pdf  of  x  is  taken  to  be 

fQ(x)  =  fL(x,d), 

which  models  linear  CFAR,  with  L  samples  in  the  threshold,  but  with  signal 
components  (unwanted  in  the  CFAR  application)  included  in  these  samples.  Only 
the  sum  of  these  signal  powers  affects  the  pdf,  and  the  average 
signal-to-noise  ratio,  per  threshold  sample,  will  he  d/L. 

In  further  generalizations,  a  is  allowed  to  fluctuate,  using  the  same 
models  as  in  Section  4  and  also  (separately  and  together)  the  new  threshold 
signal  components  are  randomized  using  the  pdf 


(1-12) 


fo<d>  -  K  H(b  - 


for  the  parameter  d.  The  average  SNR,  per  threshold  sample,  in  this  last  case 
is  Kh/L. 

A  key  feature  of  the  basic  approach  is  the  duality  between  fluctuating 

signals  and  fluctuating  thresholds.  In  Section  3  a  useful  artifice  is 

introduced  which  allows  us  to  obtain  the  detection  probability  for  a  fixed 

signal  and  threshold  fluctuating  according  to  a  pdf,  fQ(y),  from  the  solution 

of  the  dual  problem,  where  y  is  constant  and  a  fluctuates  with  the  same  pdf, 

f  (a), 
o 

The  results  are  collected  without  proof  in  Section  2,  and  the  basic 
method  is  developed  in  Section  3.  The  standard  formulas  are  obtained  in 
Section  4,  and  some  extensions  of  these  results  are  given  in  Section  5. 

Finally,  in  Section  6,  some  alternative  finite  forms  are  obtained  and 
for  ore  common  case  a  derivation  of  the  finite  fora  directly  from  the 
equivalent  infinite  series  is  given. 

n.  COLLECTED  RESULTS 

The  results  of  this  study  are  collected  here  without  proofs,  all  of  which 
will  be  found  in  the  next  two  sections.  In  all  cases,  the  assumptions  and 
conditions  are  spelled  out  in  detail,  but  detection  and  false  alarm 
probabilities  are  called  simply  and  P„A,  since  the  number  and  variety  of 
parameters  become  too  great  to  list  as  arguments  in  some  of  the  cases.  The 
following  conventions  are  used  for  the  most  coimnon  parameters: 


Parameter 

N 

L 

M 


Significance 

Number  of  "signal-plus-noise”  samples 
integrated  for  detection. 

Humber  of  "noise-only"  samples 
Integrated  to  establish  a  CFAR  threshold. 

Signal  fluctuation  parameter,  corresponding 
to  chi-squared  distribution  with  2M  degrees 
of  freedom. 
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K 


Fluctuation  parameter  for  signal  components 
in  threshold  samples,  again  chi-squared,  with 
2K  degrees  of  freedom. 

a  Signal  parameter  for  fixed  signals 

b  Signal  parameter  for  fluctuating  signals 

c  CFAR  threshold  multiplier 

d  Parameter  for  fluctuating  signal 

components  in  threshold  samples 

Three  functions  enter  these  formulas.  These  are  defined  here  and 
discussed  elsewhere. 


(1)  PN(y,a)  =  /  f„(u,a)  du  ,  (2-1) 

y 

where  f,.  is  given  by  (1.2).  P„(y,a)  is  the  probability  of  detection  for 

i»  a 

fixed,  normalized  total  signal-to-noise  ratio,  a.  It  is  essentially  the 
Marcum  O-function  for  incoherent  integration  of  K  samples. 


(2)  ss(,)  -  PK(y,0)  -  e'*  Y  k 

m=0 


(2-2) 


S,,(y)  is  the  false  alarm  probability  corresponding  to  ?„(y,a).  It  is  computed 

a  a 

recursively: 


s,„,  (y)  =  s,.(y)  +  T„(y) 


1 


N-l 


(3)  RN(y»M)  =  l  (M"J;+n)  yn  . 


n=0 


These  functions  are  also  computed  recursively: 


RN+!(y,M>  =  Vy’M)  +  Vy*M) 


VN(y,M)  -  (1  +  ~)y  VN_1(y,M)  / 


R^y.M)  =  Vq  (y,M)  -  1 


N>1 


(2-4) 


(2-5) 


and 


The  probability  density  functions  used  for  the  flu  -uation  models  are 


N-l 

%<*>  ■  (fejr  ^ 


N-l 

2 


fN(x,a)  »  e  (-)  1n_1(2/ox) 


A.  Non-Fluctuating  Signal 

In  these  formulas  the  signal  components  of  the  random  variable  being 
tested  for  signal  presence  are  constants.  The  average  SNR  per  sample  in  all 
these  cases  is 

SNR  =  a/N. 


(1)  Fixed  Threshold,  y 
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(4)  Fluctuating  Threshold,  y=cx,  conditional  pdf  of  x  is 
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f  (x | d)  ■  f  (x,d)  ;  pdf  of  d  is 

O  L. 


f0(d)  "  ¥  gK(f)* 


Condition:  K>L 


p  „  !  _  cN(l+h)N+L~K  KyL  (K-L  (__h_)J 

(l+c+ch)^-1  iio  *  1+C+Ch  ' 


*  Vl  <I*£h> 

n=0 


p  .  (l+h)L~K  ___  KyL  (K-L  ( h  * 

FA  d^h)^-1  tio  *  1+c+ch 


V  ,N+L-l+£v  ,  .  ,  .n 

x  l  (  )  (c+ch) 


B.  Fluctuating  Signal,  pdf  of  a  is 


V’  “  b  8M^  b  ^ 


Average  SNR  per  sample : 


(1)  Fixed  Threshold,  y,  P„,  given  by  (2-7) 


(2-12) 


(2-13) 
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M>N: 


PD  - 


Cl+b) 


1  M"N  U-M  . 

<  A  >  b  SN+A  (l+b> 


M<N: 


PD  “ 


(^)N-M[?1  («:* 


'>  (-b)‘  Vife) 


N-M  A  | 

£  <Vl‘>  <ife>  Vy)| 


Special  cases: 


M-N: 


PD  "  SN  (I+b) 


M-2N: 


PD  " 


ci+b)N  io  <‘)  bt  Swt 


M-l  (N>2) : 

1+h 

•  [i 


'  SN-1  (&)] 


M-2  (N>3) : 


l+b,N"2  "  1+b 


PD  *  SN-2<*>  +  <T>  • 


|a  -  Nb2)  u  -  sN_2(iTb)J 


(2-14) 


(2-15) 


(2-16) 


(2-17) 


(2-18) 
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. . . . . . . . . pip* . . . . . . . . . . . . . . . . . . . 


+ 


y 

l+b 


(2-19) 


M*2N: 


l  (")  b1  .  « 


N 


D  (l+b+c)L 


.  a 

■  V  0  “N+-£v  1+b+c 


(2-23) 


12 


ti»l: 


(1+b) 


N+L-l 


N-l 


D  bN-1(l-.b+c)L 


(He)1*"'1  £^0  '  4 


"  1  ‘  i+crrc 


M=2: 


(14,b)N^.L-1 -  jl  -  (N-2)  ^  +  (N+L-l)  C 

L+l  b 


D  bN_2(l+b+c) 


1+b 


w  4-  ii+<£+2-N)i^1  ^ 

)  1=0  ( 


(1+c) 


£+2- 


(3)  Fluctuating  threshold,  y=cx,  pdf  of  x  is 


f  (x)  =  f  (x,d). 
O 


P_,  given  by  (2-11) 
FA 


M>N: 


D  d+b+c)84-1-'1  A  *  1+b+c 


x 


N+Z-l 

I 


n=0 


N+L+  £-1 s , 1+b s  c  ,  cd  ^ 
■  L+n  K  c  }  n+l4+b+c; 
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Special  cases: 


M=N : 


„  _  cN  ^l+b)1  V  ,N+L-lx/l+b,n  e  ,  cd  _x 

P~ - IT^*  1  {  L+n  K  c  ;  n+ri+b+c; 


(Ifb+c)^-1 


n=0 


11=2K: 


p  . 

D' 


*■  V  Cl+b+c' 


£=0 


N+£~l  TIL11  J 

r  /N+L+?-1n ,l+bN  „  ,  cd  x 

*  '  L4r.  )(— >  WltS*’ 

n"0 


li*l : 


P  =  (l+b)^"1  1+c+d 

D“  b^dd-b^)1^ 


N-l  H-l  M,T  , 

c  i  x 


„  .  XN+L-1  ^  '  L+£ 

(1+c)  £=0 


cd 


,  ,lx£  o  ,  cdx  1+c+d  ,l+b+cx  ,,  r  bed 

x  *  (c)  s£+i(t+c}  ■ 6  (-bT~)  S£+1  1  a+eTcr+b+c) 


(4)  Fluctuating  threshold,  y=cx,  conditional  pdf  of  x  is 

f  (xid)  ■  fT(x.d)  ,  pdf  of  d  is 
O  Li 

Vd>  ■  ¥  8k<¥> 

Any  H,  K>L: 


(2-27) 


(2-28) 


(2-29) 
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,  cVjj.  )atL-K(14c4.h)H-|<-1-+1  x 

P  =  1 - - 

"  (1+b+c+ch) 


K;L  ,K-LV  h  *  Lr£_1  .IH-L+t-iw  ^h)n  R  (  b  ,  M)  (2-30) 
I  <  «  )(1+r+7h)  L  (  N+n  Hc+ch;  n+lk 1+b+c+ch 


*  £- 0  '  *  /Vl+C+Ch  n=0 


PpA  given  by  (2-13) 


Any  K,  M>N: 


PD' 


(1+b+c+ch) 


K 


M;N  ,M-HW  be  *  N+r_1  ^L+A-lwi+k^  r  ( - -  ,  K)  (2-31) 

x  I  ^  e  ><T+b+T>  '  '  L+n  ^  c  '  n+1 '‘1+b+c+ch 

£=0  n=0 


N-l nK-N-L+1  N-l 
c  (1+c) 


'FA 


(1-fc+ch) 

Special  Cases: 

K=L,  any  M: 


K 


rx  ,N+L-1.  ,1  »  „  /  ch  Rv  (2-32) 

X  <  L+n  }  (c>  Rn+l  l+c+cir’  K)  U  ^ 


n=0 


.  AihwW)"^ 

P  =  1 - - - - w 

D  (1+b+c+ch) 


L-l 


V*  .N+L-lx,  ^.xii  „  ,  b 

I  (  w+n  Rn+1  1+b+c+ch 


M) 


n=0 


(2-33) 
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(2-34) 


P 


FA 


1 


. N+L-l 

1.1-1  c-i-cn; 


N— 1 

I  ( 

u— 0 


N+L-l 

n 


)(c+ch)n 


I1=N,  any  K: 


P 


D 


N-l 


(1+b)  (1+b+c) 


K-N-L+l 


( 1+b+c +ch)^ 


N-l 

l 

n=0 


. N+L-l v 1+b' 
1  L+n  n  c  ; 


n 


R  ^  ( 


ch 


h+1  1+b+c+ch 


K) 

(2-35) 


PFA  given  by  (2-32) 


C.  Fluctuating  signal,  pdf  of  a  is 


f  (a)  =  i  fM(~ 
o  o  a  b 


d)  . 


Average  SNR  per  sample:  SNR  =  ^  (!l+d) 

(1)  Fixed  threshold,  y.  Pp^  given  by  (2-7) 
ii  >N : 


PD  = 


(1+b) 


M-N 


ti-N 

l 

£=0 


(li-N) 
v  l  ' 


b  P. 


( 


y 


bd 


N+£  vl+b  ’  1+b 


) 


(2-36) 


It  is  often  required  to  evaluate  detection  performance  when  the  false 
alarm  probability  is  fixed,  which  involves  the  inversion  of  the  formula  for 
PpA  in  order  to  determine  the  threshold.  The  Newton-Raphson  iterative 
method(^)  is  useful  for  this  purpose,  and  we  discuss  its  application  to  the 
two  most  common  false  alarm  formulas,  (2-7)  and  (2-9).  The  other  cases  listed 
above  involve  CFAR  problems  with  signal  components  in  the  thresho  d  samples, 
but  these  components  are  usually  not  anticipated,  and  the  threshoi.  will  have 
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been  found  from  (2-9).  Formulas  like  (2-11)  and  (2-13)  then  show  the  effect 
of  the  signal  components  on  PpA  after  the  threshold  has  been  fixed. 

Beginning  with  (2-7),  we  want  to  solve  the  equation, 

$(y)  2  SN(y)  -  ppA  -  0, 

for  y,  where  PpA  is  the  assigned  probability  of  false  alarm.  The 
Newton-Raphson  iteration  is 


y 


n+1 


0<vn), 


where 


G(v) 


$(y) 

$'(y) 


and  yc  is  a  suitable  initial  value.  A  good  plan  is  to  start  the  iteration 
at  a  point  where  the  derivative,  $’(y),  is  large,  to  avoid  wild  oscillation  of 
the  sequence,  y^.  This  can  be  assured  by  choosing  vq  to  be  the  solution  of 

♦-(yQ>  =  o. 


which  will  be  unique  in  the  applications  made  here. 
From  (2-2)  we  have 


$’(y)  “  SN(y)  =  -e 


-y 


N— 1 

m 

N-2 

m 

V 

Ve'" 

m! 

y 

y_. 

a 

ii  f 

o 

m=0 

m! 
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In  C 2—3)  we  nave  written,  in  effect, 


N-l 

s„(y)  -  I 

n=0 


where 


» 


hence 


SN(y)  =  "Vl(y)* 


As  Sy(y)  is  computed  recursively,  from  (2-3),  for  use  in  the  iteration,  the 
final  term,  T^_^(y),  can  be  recovered  for  the  evaluation  of  $'(y).  Then 


G(y) 


y  + 


sN(y)  -  pfa 

Wy) 


> 


and  iteration  is  continued  until  the  correction  term,  G(y)-y,  reaches  some 
preassigned  small  value.  We  also  find 


r(y ) 


\ 1 


N-l 

[y  -(N-l)y 


> 


and  so 


y  =  N—l 
o 


is  the  desired  starting  value. 


A  very  similar  procedure  can  be  used  for  (2-9)  with 


1  r  ,  N+L-l  s  _.£  „ 

4>(y)  o^nN+L-1  -  *  £  5  *  "  ‘FA* 

(1+y)  £=0 


(we  have  retained  the  symbol  y,  instead  of  the  c  of  ( 2—9) ) . 
This  time. 


\  N+L-l  r  ,m-l,  .  £ 

»  ~  ~  ■■  .IM.  .1  '  l  )  y 

(1+y)  £=0 


N— 1 

,  1  V  (N+L-l)?  £-1 

(l+y)^1’1  £ 1  (N+L-£~x)  !(£-!) !  > 


-^lY  C*^1)  y'  -  Y  (■£*><,«  ♦  yS 

(l+y)‘  L  <  £=0  £  £=1  1  1 


Using  the  reduction  formula 


, N+L-l ,  =  -NbL-2.  ,N+L-2. 

^  £  ’  '  £  }  K  £-1  1  ' 


the  quantity  in  curly  brackets  becomes 


. . .  ii[  i|(iip  i«  win  ii‘  iiimii  +*  iiiilll|lll||lll|||i|||iilHIIIIIIIIIII”ill|l^lllii 


N-l 


UiT  O 


N-l 


C  *  r.  ,HTL-4,  X 

1  (  c  )  y  -  >  ^  ,  )  y 

£=0  X  £=1  *  1 


»»*]  *» 
N-i }  y 


Therefore 


♦*(y)  =  - 


N+L— 1  , N+L— 2 »  N— 1 

■<  )  y 


(l+y) 


N+L 


L  .  1  .N+L— 1-  N-l 

1+y  Cl-K')^"1  !'_1  7  ' 


and  if  we  define 


n  fy)  =  _ 1 _  "y1  „* 

V  4-0  ‘  J 


N-l 


=  I  Q»(y)» 
£=0 


then 


q£(y> 


1  ,N+L-1»  £ 

<  *  >  y 


CI-^-1  * 


and 


*'(y) 


-1%  Vl^- 


Thus  the  iterative  algorithm  is  based  upon 
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G(y>  = 


wN(y)  -  p 


FA 


5 


and  ,  (y)  is  recovered 
it— i 


fron  the  iterative  coaputation  of  W  (y): 


w  (y) 

XI 


Vl<*>  +  Wy) 


,(y)  =  -  l)y  Q£_2(y) 


^(y)  -  Q0(y) 


(i-*y> 


N+L-i 


Finally,  we  compute 


/ U+L— 1 ,  (K-l)yK~2(H7)-(lfrt.)y"  1 
N“1 


and  therefore 


v 

'O 


L+i 


provides  a  suitable  initial  value, 
in,  BASIC  FOHHULATIOS 

Tna  whole  analysis  here  is  based  on  Schlafli’s  integral 

(-Q\ 

't  for  the  Bessel  functions.  We  use  it  in  the  form 
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In(x) 


n  i 
v2;  2tt1 


2 

x 

4t  dt 


t 


n+1 


(3-1) 


The  contour  of  integration  is  a  small  circle,  enclosing  the  origin  in  a 
positive  sense.  With  this  representation,  the  pdf  of  u  (for  fixed  a). 


N-l 


fN(u,a)  =  e  u  a(J)  2  1^(2®, 


(3-2) 


can  be  written 


t  + 


au 


c  ,  s  -u-a  N-l  1  r  t  dt 

V“'a)  =  e  “  sr  k|.ee  7 


Replacing  t  by  ut  in  the  integral,  we  obtain 


fN(u,a)  =  e 


-u-a  1 
2iri 


ut  +  — 

t  dt 


|tl=e 


N  * 


(3-3) 


This  representation  is  more  convenient  than  the  Fourier  transform 
expression.,  which  can  be  obtained  from  it  by  the.  change  of  variable  t=l-iX. 

We  carry  out  this  transformation,  since  we  need  the  characteristic  function  of 
f^  later  or..  After  substitution,  we  have 


fN(u,a)  =  e 


1 

2ni 


-iuX  + 


/  e 
C 


a 

1-iX 


(-idX) 

(l"iX)N 
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where  C  is  a  small  circle  enclosing  the  point  X  *  -i.  Tills  contour  can  be 
expanded,  without  changing  the  integral,  until  it  runs  along  the  real  X-axis 
from  +A  to  -A,  and  then  back  along  an  arc  of  radius  A,  in  the  lower  half- 
plane,  until  it  closes  at  X  *  +A  again.  Because  of  the  factor  exp(-iuX)  in 
the  integrand,  the  integral  along  the  arc  vanishes  in  the  limit  A+«,  for  any 
N=l,2,....  We  then  reverse  the  path  along  the  real  axis  to  obtain 


fN(u,a)  =  e 


00 


/ 


-iuX 

e 


a 

1-iX 

e 

(l-iX)N 


dX 

2tt 


The  inverse  transform  gives  the  characteristic  function 


«  .  „„  1-iX 

/  e  f„(u,a)du  =  e~n  - -  =  fc/X.a)  • 

0  N  (1-iXr 


(3-/.) 


Recalling  the  definition  of  a  in  equation  (1-3),  we  can  write  <j>N(X,a)  as  a 
product: 


¥x) 


N 

7T 

n=l 


-a 


n 


_ n 

1-iX 

e 

1-iX 


2  2 

where  an  =•  Is^l  /2a  ,  which  corresponds  to  the  definition  of  the  random 
variable,  u,  as  a  sum.  Each  factor  here  is  the  characteristic  function  of  a 
non-central  chi-squared  random  variable,  with  two  degrees  of  freedom,  and  u  is 
non-central  chi-squared  with  2N  degrees  of  freedom. 

From  our  basic  integral  representation  for  f^(u,a)  we  obtain  the  "no 
signal"  special  case, 
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This  set  of  pdf’s  will  be  used  to  model  fluctuating  thresholds  and  signal 
parameters  later  on,  and  we  introduce  the  notation 


SN<U>  =  Vu>°)  • 

Evaluating  the  contour  integral  for  f^(u,0),  we  get 

/  X  _  1  N-l  -u 

8N(u)  “  (N-l)!  u  e  * 


(3-5) 


which  is,  of  course,  the  chi-squared  pdf,  for  2H  degrees  of  freedom. 
Obviously,  Eu«N,  if  g^(u)  is  the  pdf  of  u. 

The  "probability  of  detection"  corresponding  to  the  pdf  fN(u,a)  is 
denoted  ?^(y,a)  and  it  is  defined  by  the  equation 


PN(y»a)  =  /  fN(u,a)du  . 

y 


(3-6) 


Our  notation  differs  from  the  standard  one  (4)  in  that  we  retain  the  same  order 
of  the  variables  as  in  the  pdf,  and  the  parameter  a  is  the  "total  SNR,"  as 
defined  Ln  Section  1.  The  advantages  of  this  choice  seem  worth  the  cost  of 
departing  from  common  usage.  In  terms  of  P„,  Marcum's:  Q-function  is 


To  obtain  an  integral  representation  for  P^,  we  substitute  the 
representation  (3-3)  in  definition  (3-6)  and  reverse  the  order  of 
integration: 


PN(y,a>  -  e'a 


(  t 

It  -  8 


t  -(l-t)u  dt 

J  e  du  17 


Since  Re{l-t}  >  0  on  the  t-contour,  the  u-integral  is  uniformly  convergent, 
thus  validating  the  interchange,  and  we  obtain 


?N(y,a)  =  e 


-y-a  1 


yt  +  t  dt 


1 1 1  =*e 


tN(l-t) 


(3-7) 


This  is  the  fundamental  representation  for  P^,  used  throughout  this  study. 

A  second  representation  is  obtained  by  making  the  change  of  variable,  t+l/t, 
in  (3-7): 


pN(y»a)  -  e 


-y-a  1 
“  6  2iri 


I  1 1  =1  /  G 


*t+£  t^dt 

5  t-1 


The  contour  now  is  a  large  circle,  including  the  simple  pole  at  t=l,  as  well 
as  the  essential  singularity  at  the  origin.  The  residue  of  the  integral  at 
the  simple  pole  is  exp(a+y),  and  hence  we  can  write 
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„  ,  -y-a  1  ,  at  +  t  tN_1dt 

P^Cy.a)  -  1  -  e  Til  ,  e  "Ft 

I  t|=e 


(3-8) 


where  the  contour  is  again  a  small  circle  enclosing  the  origin. 

Both  representations  are  useful  and  we  note  that  the  integral  in  (3-8)  is 
just  like  the  one  in  (3-7),  except  that  the  variables  y  and  a  are 
interchanged,  and  the  factor  t  appears  with  a  non-negative  exponent.  In  other 
words,  we  have  the  evaluations 


-y-a  1 


yt  +  t  t"ndt 


?n(y,a)  ;  n>0 


1-Pln(a,y)  ;  n<0 


It  proves  to  be  exceedingly  useful  to  define  PN(y,a)  for  all  integral  N  by 
means  of  representation  (3-7),  with  the  understanding  that 


P_»(y»a)  si-  P  (a,y). 


(3-9) 


This  artifice  allows  us  to  exploit  the  duality  between  threshold  and 
signal  parameter  In  the  following  way.  Suppose  the  signal  parameter  is 
randomized,  according  to  some  pdf,  fQ(a)*  With  fixed  threshold,  y,  the 
probability  of  detection  will  be 


FM(y)  =  I  PM(y,a)f  (a)da. 


(3-10) 
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Now  suppose  the  right  side  is  evaluated  for  positive  and  negative  N,  using 
(3-7)  as  the  definition  of  P^(y,a).  Then  the  dual  problem,  in  which  a  is 
fixed  and  y  is  random  (with  the  same  pdf)  has  detection  probability 


00  CO 

/  PM(y,a)f  (y)dy  =  /  U  -  P  „(a,y)]  f  (y)  dy 
0  N  °  o  1  N  0 


=  1  -  Fj_N(a)  •  (3-11) 

This  technique  saves  the  needless  duplication  of  calculations,  which  would 
otherwise  use  (3-7)  for  one  set  of  problems,  and  (3-8),  in  a  parallel  way, 
with  the  dual  problems. 

The  source  of  this  duality  is,  of  course,  the  basic  pdf.  The  fact  |hat 
ln(x)  -  I_n(x)  can  be  proved  from  (3-1)  by  making  the  substitution  t 
in  the  integral  and  hence  from  (3-2), 


Wu’a)  =  e 


-u-a 


-a-u 

e 


N 

2  F^C^au) 


Thus  we  could  define 


fl-N(a,u)  =  fW4-i(u’a)* 


m' 


(3-12) 
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and  this  definition  is  consistent  with  (3-6),  extended  to  negative  values  of 
N.  Integration  of  both  sides  of  (3-12)  over  a  yields 


/  f1_N(a»u>da  3  Pi-N(b’u)  =  1  "  V'u»b) 

b 


°  /  fN+i(u,a)da’ 

b 

or 

00 

PN(y,a)  =  1  -  /  fN+1(y,a')da',  (3-13) 

a 


an  interesting  result  which  also  follows  directly  from  (3-3). 

The  probability  of  false  alarm  for  fixed  threshold  will  be  denoted 

SN(y)  =  Vy’0)* 


From  (3-7)  we  have  the  integral  representation 


SN(y)  3  ^  2 ¥1 


[  eyt  —gi¬ 
lt  I  =e  tN(l- 


t) 


(3-14) 


which  will  often  be  used  to  ivaluate  contour  integrals  in  this  analysis.  For 
N>1,  we  have 


00 
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after  expanding  (1-t)  \  and  then 


N-l  n 

SM(y)  =  e  ^  £  ^~r  • 

N  J  Ln  m! 

m=0 


(3-15) 


The  series  terminates,  since  there  is  no  pole  at  t=0  when  m>N.  Of  course. 


sN-(y)  "  I  SN(u)du  , 

y 


which  again  yields  (3-15)  after  substitution  of  (3-5)  and  repeated  integration 
by  parts. 

From  (3-14)  we  have 


S  (x)  ®  0  ,  n<0  . 

n 


(3-16) 


This  is  completely  consistent  with  (3-9),  since 


S_N(x)  3  P_N<X*°)  3  1  ~  Pn+1^°’x^’ 


and 


V°»x> 


3  I 


fN(u,x)du  =  1, 


by  normalization. 

We  proceed  now  with  the  general  problem  of  detection  with  a  fixed 
threshold  and  random  signal  parameter.  As  above  we  let  fQ(a)  be  the  pdf  of 
a,  and  we  define  the  characteristic  function 
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(3-17) 


j>  (X)  =  /  e^af  (a)  da 

0  0  0 


We  evaluate  the  detection  probability  by  substituting  the  integral 
representation,  (3-7)  into  Eq.  (3-10): 


/  £o<a) 


-y-a  1 
s  2iri 


I  t{=e 


yt  +  f  dt 


tN(l-t) 


da 


We  want  to  reverse  the  order  of  integration,  which  will  involve  the  evaluation 
of 


-(1  - 


I)a 


f  (a)  da  . 
o 


For  this  integral  to  converge  uniformly,  it  is  necessary  that  Re(l-  t)  >  0, 
which  will  be  satisfied  if  we  can  arrange  to  have  ! 1 1  >1 .  This  is  accomplished 
by  expanding  the  original  contour.  It  I  =  e,  until  it  becomes  a  circle  with 
radius  larger  than  unity,  and  compensating  for  the  for  the  effect  of  the  pole 
(at  t=l)  now  included.  In  other  words  we  write 


-y-a  1 
2ui 


|t|=e 


yt  +  - 


dt 


t  (1-t) 
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=  1  +  e 


-y-a 


1 

2iri 


ltl>l 


yt+f 


dt 


tN(l-t) 


and  Chen  reverse  the  order  of  integration.  The  result  is 


FN(y)  =  1  +  e 


2ni 


[  t|  >i 


00  -(1-  ■— )a 
/  e  f  (a)da 

is  o 


dt 


N.,  , 

t  <l-t) 


The  integral  over  a  yields  $  [i(l-  ■£•) } ,  with  uniform  convergence  to  justify 
the  interchange. 

Now 


UD(X)I  <  1 

whenever  Im(A)  >  0,  and  this  upper  half-plane  maps  into  the  exterior  (and 
boundary)  of  the  circle. 


It 


1 

2  » 


in  the  t-plane.  It  will  be  true  in  the  cases  considered  here  that  4>0(A)  is 
analytic  when  Im  A  >  -  e,  for  a  suitable  positive  e,  so  that  <j>o(i(l  -  ^-) ] 
will  remain  analytic  for  a  small  but  finite  distance  inside  the  circle 

It  --^1  “  j.  This  will  allow  the  t-contour  to  be  shrunk  just  inside  the  unit 

circle.  Since  <f>Q (0 )  =  1,  the  pole  at  t=l  is  simple,  and  hence 
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dt 


F  (y)  =  1  +  e'y  f  eyt  0>  [i(l  -  - 

ivl  | t I >1  0  C  t  (1-t) 


or 


V*> 


-y  1  r  yt  ,  rw,  lN,  dt 

*  e  ■=— r  j  e  $  [i(l  -  -)]  — - 

2tti  if,  o  t  N,,  . 

I t I =l-e  t  (1-t) 


(3-18) 


This  is  the  basic  expression  which  will  be  used  throughout  this  report. 

IV.  DERIVATION  OF  THE  STANDARD  FORMULAS 

We  begin  by  randomizing  the  signal  parameter,  using  a  general  even-order 
chi-squared  pdf  for  a.  In  particular,  we  postulate 


fo'a>  -  I  %<!> 


(4-1) 


where  g  (x)  is  the  pdf  defined  by  Eq.  (3-5). 

il 

then 


The  characteristic  function  is 


,  /.\  f  1^3  ,aN  da 
*o(X)  =  J  e  &,(-) 

0 


^i'b'  b 


=  J  elbXx  g,,  (x)dx 
0  11 


-  a  -  ibxr1  , 


a  special  case  of  (3-4),  namely  <^(bA,0).  With  this  model, 

Eu  =  Ea  +  N  =  Mb  +  N; 

and  if  tl  =  kN,  the  result  is  the  same  as  assigning  to  each  component  signal 
parameter,  a  ,  the  pdf 


f(V  -  £ 


y 


with  each  component  independent  of  the  others.  In  this  case,  Ea  =  kb.  The 

n 

Swerling  cases  1,  2,  3  and  4  correspond  to  the  M-values  1,  N,  2  and  2N, 
respectively. 

We  note  that 

•otl(1~  I5 3  =  [  (l+bH-b1 

and  abstitute  in  Eq.  (3-18).  The  resulting  expression  for  the  detection 
probability  is 


-y  1 


2iri 


lt|=l-< 


,yt 


M-N 


dt 


f(l+b)t-b]h(l-t) 


V 


4-2) 


A  basic  difference  appears  now,  depending  on  whether  11  w*  or  !1<N.  We  suppose, 
first,  that  it>N,  so  the  only  pole  in  (4-2)  is  at  t=b/(l+b).  We  shrink  the 
contour  to  a  small  circle  enclosing  this  pole  and  then  make  the  change  of 
variable 
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. . . . . . . . 


(4-5) 


the  desired  expression. 

The.  false  alarm  probability  is,  of  course- 


PFA  -  V*> 


as  it  must  be  for  any  fixed-threshold  model.  Note  that  our  pdf,  fQ{a), 
approaches  6(a),  as  h*0,  and  ^Q(^)  +  1*  One  can  also  put 


b  =  a  /It 
o 


and  let  M+«.  In  this  case  f  (a)  aporoaches  6(a-a  ),  and  6  (X)  exp(ia  X). 

o'  o  o  o 

while  Pp  is  transformed  into  a  standard,  infinite-series  expression  for 
PN(y,ao)  (this  is  Fehlner's  series,  which  is  discussed  in  Section  6). 

For  fi=N  we  have  the  familiar  Swerling  case  2  result 


PD  “  SN^l+b^  > 


so  favored  by  radar  analysts.  When  H=2!i,  we  obtain  the  Swerling  case  4 
expression 


pn  =  jj  I  Ab1  C*-2») 

D  (1-h>)n  £= o  1  1+0 


Whesi  M<N  we  write  (4-2)  as 


P  =e"y  -i 
D  2x 


/  ^ 


I  ti=l-e  ul+b)t-b)  t  (1-t) 


fa-fi'i 
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and  proceed  differently.  The  integrand  now  has  two  poles,  and  we  write 


p  =  P  +  P 
D  rl  2 


(4-7) 


whece  the  integrands  of  P^  and  P^  are  the  same  as  PQ,  but  the  contour  for  P^ 
is  a  small  circle  enclosing  the  pole  at  t=b/(l+b),  while  for  it  is  the 


circle  |t|=e. 


The  evaluation  of  P^  is  the  same  as  the  case  just  analyzed,  which  used 


transformation  (4-3),  up  to  the  point  where 


P1 


_  _L_ 

l+b  1 


ds 


ttt  s  ■  •  -(N-ti) 

/  l+b  ,s+bs. 

2 it:.  -  |  6  4+bJ  It..  N 

|s|=e  s  (1-s) 


This  time  the  expansion  is  an  infinite  series: 


-(N-it)  ,  «> 

/■s+bx  /-l+b,,  r  ,tHi-l+JL  ,  .  .1 

(l+b)  ("TT>  L  (  #  )(-s/b)  , 


fc=0 


which  converges  since  |s|=t'  and  e  can  be  chosen  smaller  ‘.nan  b.  When  this 
series  is  substituted  and  the  integration  performed  termwise,  only  the  terms 
with  £<J1  contribute,  since  beyond  that  there  are  no  poles.  Theiefore, 


n-m  M-l 

p  -  A+-V  V  fN-M-l+£u_  ,  Jl  .  y  . 

1  '  zl0  *  K  7  Wl?^  ’ 


(4-8) 


where  (3-14)  has  again  been  employed. 


i  . 
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can  be  handled  in  a  similar  way,  without  changing  the  variable  of 
integration  from  that  used  in  (4-6)  but  expanding 


KUMt-bl""  -  (-  i)"  i  («-l+*)(«t)‘. 
J  £=0 


Again,  since  |t|=e  ,  convergence  can  be  assured  and  we  get 


P  -  (-l/b)"  ?  /  «>'  _ «£ _ 

2  \  t-ruJ  2iri  .  <  N-m-j 

£=  0  *  D  I t|=e  t  X(l-t 


This  time  the  sum  stops  at  £=N-il-l,  and  the  integrals  are  evaluated  using 
(3-14).  Altogether,  after  reordering  the  suras  to  get  an  increasing  sequence 
of  S-f unctions, 


M  1  4-k  C 

PD  ■  P1  +  p2  -  (-1/b>  <TT> 


f  y  (b'1-|1k- 
I A  !H1-1 


b>* 


N— 11  „  ,  B  ,  £  1 


(4-9) 


The  usual  application  of  these  formulas  will  be  to  Swerling  cases  1  and 
3,  where  11=1  and  11=2  respectively.  For  these  or  any  small  values  of  11, 
another  expression  is  more  convenient.  It  is  obtained  from  (4-6)  by  writing 


(l+b)t-b 


and  put 


the  complete  expression  is 


wy> 


(l-x)M 


01-1) ! 


(*— ) 
vdx' 


11-1 


H 

(d- 


(x-l)y 


x)x 


N-M 


[1  "  SN-M(xy)1 


(4-i: 


where 


x 


b 

1+b 


The  identity 


df  6X7 


y  exy [1  -  S^Cxy)]  , 


which  is  easily  verified,  makes  the  evaluation  of  the  derivatives  in  (4-11) 
relatively  straight-forward. 

For  Swerling  cases  1  and  3  we  have: 


N-l - 4— 

PD(tl-l)  =  S^_1  (y )  +  (i±t)  e  1+b 


[1  SN-ld+b^  ’ 


(4-1? 


and 


Pd(M-2) 


W>->  +  (ir)N"2  ■ ' lib  |a  -  [i  - 


+  i%  t1-^ 


K><Sb»j 


(4-13 
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By  writing 


S  (JS-0 
N-ri+b' 


Vife>  - e 


by 

l+b 


N-l 


(N-l) ! 


(ZL -) 
4+b' 


and  substituting  into  (4-12),  it  follows  that  this  equation  can  also  be 
written 


yy> 


<¥> 


N-l 


y 

l+b 


[i  -  S, 


4+b'J 


(4-12) 


In  this  form,  (4-12)  is  also  valid  for  N*l,  in  which  situation  Swerling  cases 
1  and  2  coincide.  By  a  similar  manipulation,  (4-13)  can  be  written 


SN-l(y)  +  e  7  (l+b)(N--2) ! 


N-l 


+  (~) 


N-2 


l+b 


(1-^ 


lfe> 


x  t  1  '  SN-l(I+b)1  (4'13) 

which  happens  also  to  be  valid  for  N=2,  where  Swerling  cases  2  and  3  coincide 

We  turn  now  to  the  dual  problem  of  fluctuating  threshold,  as  in  linear 
CFAR  detection.  The  signal  parameter,  a,  is  fixed,  and  the  threshold  is 
written 

y  =  cx, 
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where  the  pdf  of  x  is  taken  to  be  g^(x).  Thus  x  corresponds  to  the  sum  of  L 
random  variables,  each  with  the  exponential  pdf 


g1(x)  =  e  x  . 


In  other  words,  x  is  the  result  of  Incoherent  summation  of  L  complex  samples 
of  noise  alone,  the  noise  samples  having  the  same  variance  as  the  noise 
components  of  the  samples,  z^,  which  contribute  to  u.  We  can  call  these 
noise  samples  x^,  and  then 


L  x1+--*+xL 


y  =  c  l  x  =  cL 


The  average  of  the  x  is  the  noise  level  estimate  in  CFAR  problems,  and  the 

X 

corresponding  threshold  multiplier  is  cL.  With  this  model,  the  pdf  of  the 
threshold  variable,  y,  is 


fo(y>  *  I  8l{c)  • 


(4-14) 


and  the  required  detection  probability  is 


PD  =  /  Vy,a)fo(y)  dy* 

0 


(4-15) 


According  to  Eq.  (3-11),  this  is  obtainable  from  the  solution  of  the  dual 
problem: 


42 


pD  "  1  -  I  ^(a.y)  fD(y)dy 
0 


-  1  -  J  F1-Hte,7>  gL(c)  ^ 
0 


This,  of  course,  is  the  problem  we  have  just  solved  and  we  have  only  to 
make  the  changes  of  variable 


N+l-N 

M+L 

y+a 

b-*-c 

in,  for  example,  Eq.  (4-2).  We  note  immediately  that  the  case  previously 
called  '1KN"  cannot  occur,  since  the  old  exponent  M-N  becomes  N+L-l  in  the 
dual  problem.  Thus  (4-4)  provides  the  desired  solution: 


1  - 


(1+c) 


N+L-l 


N+L-l 

l 

£»  0 


( ^"L-l  \  Z  g  f  a  % 

^  Z  ;  Jt+l-N  4+cy  * 


Since  ^^(x)  vanishes  unless  m>0,  only  the  terms  in  this  sum  having  £>K 
survive,  and  our  result  can  be  written  in  the  form 


1  - 


U+c) 


N+L-l 


L-l 

l 

Z=0 


.N+L-l.  Z  _  ,  a 

(  Ml  >  c  St+1(H?) 


(4-16) 
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It  may  be  verified  that  this  result  is  also  obtained  by  a  direct 
calculation,  starting  from  (4-15)  and  using  representation  (3-8)  for  P  .  This 
calculation  closely  parallels  the  derivations  of  Eqs  (3-18)  and  (4-2). 

The  false  alarm  probablity  is  obtained  by  putting  a=0  in  (4-16)  and 
noting  that  S  (0)=1: 


P 

"FA 


>J4-T  -1 

u+cr 


N-l 


o=0 


,m-L-i.  z 
(  £  >  C  » 


(4-17) 


and  (4-16)  can  also  be  written  In  the  form 


v 

‘D 


c 


(1+c) 


N+L-l 


V  w-l-i,  _  s 

0’n  i  N+Z  '  -  11  b 

1=0 


£+l'‘l+c'J 


(4-18) 


The  probability  cf  detection  when  both  threshold  and  signal  parameter  are 
random  variables  can,  of  course,  be  approached  from  two  directions:  we  can 
start  with  random  a  and  fixed  y,  then  averaging  over  the  assumed  pdf  of  y,  or 
proceed  in  the  reverse  order.  It  turns  out  that  both  methods  yield  useful 
formulas,  not  easily  transformed  into  one  other.  We  start  with  (4-16), 
averaging  over  the  assumed  pdf  of  a: 


fo(a)  b  8mV 


(4-19) 


To  carry  this  out  we  need  the  integral 


WIIII|li|llllllllll!ll|iilllllllMI 

.  ; 


/  S£+1(Ax)  &a(x)  dx 


£  .m  ”  . 

I  h  f  6  ^  %<X>  dx 

m=0  Q 


f  xm  1  r  -(1+A)x  M-l+m  , 
l  mf  OPTTT  1  e  X  dx 

m=°  o 


-V  l  • 

(1+Aj*  m=0  ®  1  X 


(4-20) 


This  expression  recurs  in  several  of  our  formulas,  and  it  should  be  noted  that 
index  £  occurs  only  as  the  upper  limit  to  the  sum.  This  fact  will  facilitate 
the  recursive  calculations  of  double  sums,  like  (4-21)  below. 


When  (4-16)  is  averaged  over  (4-19),  we  make  the  variable  change, 
a=bx,  to  obtain 


PD  ~  1  (1+c)m-L-l  (  N+£  5  c  f  S£+l(l+c)  dx’ 


and  now  (4-20)  provides  the  desired  result 


„  ,  c"  ,  1-te  >"  V  ,S!->L-1,  i  V  b  ,, 

po  •  1  “  -Tiffixi  cT+b^>  s  <  mi  )c  u  *  <4~21> 

(1+c)  £=0  m=0 


Equation  (4-21)  is  valid  for  all  values  of  It,  and  it  is  particularly 
convenient  for  Swerling  cases  1  and  3,  where  M=1  and  2,  respectively,  Th 
falsa  alarm  probability  is  still  given  by  (4-17). 

Eor  11=1,  the  ra-sum  Js  simply 


n-0 


1+b+c 


m 

)  = 


1+b+c 

1+c 


£+1. 


1  -  ( 


1+b+c 


) 


and  (4-21)  becomes 


PD  =  1- 


( 1+c  ) 


N+L-l 


L-l  xu,  ,  n  (  .  £+1  » 

r-  N+L-l .  £  t  ,  ,  b  .  f 

*;0  '  l  [ 


1  - 


(He) 


N+L-l 


N+L-l 

) 

i=K 


f N+L-l. 


l  *  (--: 


i+b+c 


£+l-N 

) 


The  £-sum  is  now  rewritten  as  the  difference  between  a  full  sum  (£  runnir 
from  zero  tD  N+L-l)  and  a  partial  sum,  where  £  runs  from  zero  to  N-l.  Th 
final  result  is  (!1=1): 


,l+bN 


N-l 


N-l 


V  fW--l,  £  ji  ,1+bj 


N-l- 


(. 


1+b 


1 


This  formula  is  convenient  because  it  has  N  terms,  rather  than  L,  and  N  is 
usually  the  smaller  number.  It  can  also  be  obtained  from  (4-12)  by  averaging 
over  y.  When  N=l,  (4—22)  becomes  simply 


,  1+b 
4+b+c 


(4-23) 


In  Swerling  case  3  01=2),  the  m-sum  in  (4-21)  is 

£ 

£  (nri-Du31  =  (1-u)  -  (£+2)u*+*  +  (£+l)u^+^  l  , 

m=0  1  ’ 

where 

_  b 

U  ~  l+b-rc 

A  calculation  which  parallels  the  derivation  of  (4-22)  yields  the  result 
(M=2): 


(4-24) 


When  N=2,  this  formula  reduces  to 


u  +  a+1) 


(4-25) 


P 

~D 


=  (■ 


1+b+c 


) 


When!* >N,  another  formula  is  more  convenient  than  (4-21).  It  is  obtained 
by  starting  with  (4-4)  and  randomizing  y,  thus  reversing  the  order  of 
randomization  that  led  to  (4-21).  The  pdf  of  y  will  be 


and  we  get 


P„  = 


(1+b)1 


II  II 

l 

£=0 


V  „  ) 


$Nf£  ^1+b^  gL(c)  d(c) 


We  change  variables  (y=cx)  and  use  (4-20)  tc  find 


PD  =  (‘ 


l+'D 


1+b+c 


) 


(1+b)1 


M-N 

l 

£=0 


< ; )  b 


£ 

3=0 


(4-26) 


The  false  alarm  probability  is  still  given  by  (4-17),  and  the  fact  that 
(4—26)  reduces  to  (4-17)  when  b=0  rests  upon  the  useful  binomial  Identify 


l 


X-l-frnu  m 

(  _  )  y 


-  (l-y)3 


(4-27) 


which  is  proved  in  Appendix  1 


When  M=N  (Swerling  2),  (4-26)  becomes 


tJ.  i-  K-l  ,  a 

P  «  C~~)  V  (L“^  )(_§_) 

D  U+b+c'  v  m  Ml+b+c'' 

m=0 


(4-28) 


For  N=I,  Swerling  cases  1  and  2  coincide,  and  (4-28)  becomes  (4-23)  in  this 
case.  Likewise,  when  N=2,  cases  2  and  3  coincide,  and  (4-28)  then  reduces 
(4-25),  as  it  should. 

For  Swerling  case  4,  M=2N,  (4-26)  becomes 


P  =  (■ 


_l+b_, 

i+b+cJ 


(1+b) 


N 


l  Qb  l  <  L  ><T& 

£=0  *  m=0  n  1+b4c 


) 


(4-2S) 


V.  GENERALIZATION  OF  THE  FLUCTUATION  MODELS 

In  Section  4,  when  a  parameter,  say  y,  was  randomized,  a  pdf  of  the  type 


was  used,  introducing  one  new  parameter.  In  this  Section  we  generalize  by 
using,  instead,  s  "signal-plus-noise"  type  of  pdf,  such  as 


This  generalizes  the  chi -squared  pdf  to  a  non-central  chi-squared,  and  adds 


another  parameter.  Tne  model  is  equivalent  to  the  postulate  y  =  cx,  where  x 
has  the  pdf  f^(x,d)  as  defined  by  Eq.  3-2. 

We  begin  by  randomizing  the  signal  parameter,  a.  We  put  a=bx,  where  x 
has  pdf  fM(x,d).  In  other  words 


f 

o 


b  fM(b  ’  d') 


(5-1) 


and  the  characteristic  function  is 


(X)  =  E  e 
o 


iAbx 


d 

1-ibX 

■d  e _ 

M  ’ 

(1-ibX)" 


according  to  (3-4).  To  apply  (3-18),  we  rpquire 


M 


(M1^1  t^  “  6  ^ (l+b)t-b ^  exp  ^  » 


'  (l+b)t-b 


and  then 


PD  "  e 


-y-d  1 


/  e^*"  exp  [ 


2irx 


td 


(l+b)t-b 


M“N 

t  at 


,M„ 


[(l+b)t-b]“(l-t)  .  (5-2 


It-Ui 


We  confine  ourselves  to  the  case  M>N  and,  following  the  previous  analysis, 
make  the  change  of  variable  (4-3).  The  new  feature  is  the  exponential  factor, 
which  becomes 


exp  [ 


td  ,  ,d  s+b.  /  d  .  ,1  bd  N 

n+bOFb1  =  PXp(s  l+b>  =  exp(T+b)  ex^e  l+b*  • 


With  the  definitions 


d  - 


bd 

1+b  ’ 
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y 

l+b 


we  have 


-y-d 


1 

2iri 


e  W 


ds 


M,..  V 

s  (1-s) 


il-N 

As  before,  we  expand  the  factor  (s+b)  ,  and,  this  time,  evaluate  the 

integrals  using  (3-7).  The  result  is 


(l+b)' 


M-N 


M-N 

l 

&=0 


(**~^)  b^  p  (  y 

w-ri+b  ’ 


bd  . 
l+b 


(5-3) 


a  direct  generalization  of  (4-4),  to  which  it  reduces  when  d-K).  When  b=0, 
only  the  term  £=0  remains,  and  the  false  alarm  probability  reduces  to  the 
expected  value 
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The  analog  of  Swerling's  case  2  (M*N)  is  the  simple  expression 


P 


N'l+b  »  1+b'  * 


(5-5) 


A  result  of  this  kind  is  to  be  expected  (when  M=N),  since  each  complex  sample 

z  ,  is  now  being  modeled  in  the  form 
n 


z 

n 


s  +  w 
n  n 


» 


where 


s  =  Jh  (s '  +  w'  )  . 
n  n  n 


The  new  "signal  components,"  s',  satisfy 


As  a  consequence, 


z  =  /l+b  (s  "  +  w  ") 
n  n  n 


where 


Elw  "I2  =  2o2 
n 


and 


? 

3 


Equation  (5-5)  then  follows,  since  the  new  effective  threshold  is  y/(l+b). 

In  the  dual  problem,  the  signal  parameter,  a,  is  fixed,  while  the 
threshold  is  randomized  by  the  statement:  y  =  cx,  when  the  pdf  of  x  is 
f  (x,d).  Then 

Li 

■  §  £L(c  >  d>  <5'6) 

and  the  detection  probability  is  obtained  from  (3-11)  and  (5-3)  by  making  the 
changes  of  variable 

N+l-N 
M+L 
b+c 
y+a  . 
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Parameter  d  is  unchanged,  since  it  plays  the  same  role  in  both  problems.  We 
obtain 


(1+c) 


1 _ 

N+L-l 


N+L-l  m,t  « 

l  r+)  he1  p 


/  a  x 

£+l-N('l+c  *  1+c'  s 


(5-7) 


which  is  a  direct  generalization  of  (4-16).  Expression  (5-7)  reduces  to 
(4-16)  when  d=0,  since  the  terms  £-0  through  £=N-1  will  vanish  in  that  limit. 

Although  (5-3)  contained  only  P-functions  with  positive  subscripts,  (5-7) 
involves  negative  subscripts  as  well.  The  latter  are  evaluated  by  means  of 
(3-9),  but  first  it  is  useful  to  apply  (3-9)  to  all  the  terms  of  (5-7),  to 
obtain  the  formula 


(1+c) 


1 _ 

N+L-l 


N+L-l  .  , 

i  (NT  >  c*  ?, 


,_cd  _a_ . 
N-ri+c  »  1+c ' 


(5-8) 


Formula  (5-8)  is  the  basis  of  later  generalizations,  since  it  is  more  compact 
than  (5-7)  and  also  contains  the  parameters  a  and  d,  in  a  "natural  order.” 

For  computation,  we  split  the  terms  with  positive  and  non-positive 
subscripts,  using  (3-9)  on  the  latter: 


Pb  '  u+=)N+L-1 


N-l 

V  /N+L-l,  £  . 
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(5-9) 
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Finally,  the  first  sum  in  (5-9)  is  reversed,  and  the  index  redefined  in  the 
second  sum  to  obtain  the  result 
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(5-10) 


To  obtain  the  false  alarm  probability,  we  put  a=0  in  (5-10)  and  note  that 

Vi(0-  ’ l-  and 
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(5-11) 


which  reduces  to  (4-17)  when  d=0. 
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As  mentioned  in  Section  1,  the  present  model  describes  the  performance  of 
a  linear  CFAR  system  with  a  non-fluctuating  signal  and  with  other  signal 
components  in  the  terms  used  to  form  the  detection  threshold.  These  signal 
components  are  usually  unwanted,  and  they  can  strongly  modif*  the  performance 
of  a  detection  system.  One  obvious  effect  of  these  signal  contributions  is  to 
reduce  the  false  alarm  probability  (by  increasing  the  threshold),  which  can  be 
seen  from  (5-11).  All  the  terms  in  that  sum  are  positive,  and  [cd/(l+c)] 

is  a  decreasing  function  of  d,  since  S^(y)  is  the  false  alarm  probability  of  a 
system  with  n  samples  and  threshold  y,  which  decreases  with  increasing  y.  The 
linear  "CFAR"  loses  its  CFAR  property  with  this  change  in  the  statistical 
properties  of  the  noise.  It  is  evident  from  (5-10)  that  P^  is  reduced  by 
these  signal  components  also. 

Since  is  the  probability  that  u>cx,  we  have 


1-Pj,  =  Prob  <  c]  . 

Thus  1-Pjj  is  the  cumulative  probability  distribution  of  the  "non-central  F" 
random  variable,  p  =  u/x.  Both  numerator  and  denominator  of  this  ratio  are 
non-central  chi-squared  variables,  and  (5-10)  provides  a  basis  for  computation 
which  is  more  convenient  than  the  standard  series,  as  given  by  Springer. (10) 
When  d=0,  the  variable  x  becomes  simply  chi-squared,  and  the  resulting 
cumulative  probability  distribution  is  a  finite  sum  of  elementary  functions, 
given  by  (4-16).  It  should  be  noted  that  all  the  chi-squared  variables 
discussed  here  are  of  even  order,  and  this  accounts  for  the  relative 
simplicity  of  the  formulas. 

We  now  extend  our  results  by  randomizing  a  and  then  d,  finally  treating 
the  case  where  both  a  and  d  are  random  variables.  We  begin  with  a,  and 
postulate  for  this  parameter  the  pdf 
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f0(a) 


£  %<!> 


> 


as  before.  Both  sides  of  (5-8)  are  multiplied  by  fQ(a),  and  integration  over 
a  is  carried  out  after  the  change  of  variable:  a=(l+c)  x.  We  get 
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(5-12) 


The  integral  is  the  probability  of  detection  for  a  fixed  threshold  and 
fluctuating  signal  parameter,  and  hence  it  is  expressible  by  means  of  the 
results  of  Section  4  with  the  following  identification  of  variables: 


N  N-£ 


y 

b 


cd 

*  1+C 


(5-13) 


However,  the  sum  in  (5-12)  may  have  to  be  separated  into  two  parts,  where 
N-£5ti  in  one  part,  and  N-£<31  in  the  other,  unless  11  >N,  in  which  case  all  terms 
are  of  the  same  type.  We  carry  out  this  procedure  only  for  the  case  M=1 
(Swerling  case  1),  before  making  the  assumption  that  II  >N. 

When  M=l,  we  write  (5-12)  in  the  form 


P 


D 


(1+0 


1 _ 

N+L-l 


.  N+L-l. 
^  ' 


cZ  X(£) 
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(5-14) 


+ 


N+L-l 

l 


£=n 


N+L-l.  i 
o  )  c 


> 


and  use  (4-12)  to  evaluate  X(£),  and  (4-4)  for  Y(£).  First,  we  recall  that 
(4-12)  itself  can  be  written 


PD  = 


s„(y)  + 

N  D 


N-l  - 


y 

l+b 


tl  " 


sO1’ 


(4-12) 


which  is  valid  for  all  N>1.  When  the  variable  substitutions  (5-13)  are  made, 
we  find  from  this  version  of  (4-12)  that 


XU)  =  S„  „( 


cd 


K-t'l+c 


->  + 


N-  £-1 


cd 


1+b+c 


[1  ' 


S  (  k 
N-Tl+c 


cd 


1+b+c 


)] 


In  the  corresponding  sum  over  £,  we  write 


N-l 


I 

£=0 


,  N+L-l 


C*  X(£)  = 


T  =s-1^  xcN-i-o 

1=0  x 


to  obtain  an  increasing  sequence  of  S-functions.  Then, 


58 


(1+c) 


r  /  N+L-l .  ,  1*  £  /Cd  v 

\m-i  fn  1  h+z  nc}  )  i+ n+c- 


, 1+b+c x£  „  1+b+c  fT  _  c  ,  b  _ cd  S1 

+  (___)  e  [1  ^+1(— 


N+L-l  ... ,  . 

+ - kr  l  (  z  >  c  Y(£)  • 

(1+C)N  1  £=N 


To  obtain  Y(£),  we  use  (4-4)  with  the  substitutions  (5-13)  and  '1=1: 


..  1"N+*  1-N+£  .  V  rn  , 

Y(£)  =  (1+b+c)  I  (  m  )(l+£)  SN+m-£(l+b+C)  * 

m-u 


In  (5-14),  Y(£)  enters  only  in  terras  where  £>N,  and  also,  the  S-f unction  is 
zero  unless  m>£-N.  This  occurs  in  only  one  terra  of  the  ra-sum  (the  upper 
limit)  hence 


!-N+£ 

Y(£)  =  ^ 1+b+c 5  Sl( 1+b+c ^ 


Of  conT-se 


=  e 


and  we  obtain 


N-l  j  N-l 

(l+c)^-1  I  Jo 


,N+L-1. ,1*£  /Cd_\ 

K  L+£  'V  b £+1 ' 1+c ' 


cd 

1+b+c 


£f0  (  L+£  H  be  )  U  S£+ri+c  1+b+c ' 


~  1+b+c  ^  ,N+L-1W  be  \  £+H 

+  e  Jo  (  »«  >(W  {• 


(5-16) 


The  terms  not  involving  S-functions  in  (5-16)  can  be  combined,  as  follows 


Nr2  ,N+L-1W  1+b+c/,  Lr1  ,N+L-1W  be  /+1 
L  (  L+£  H~bT~ )  +  i  (  M+9 


N+£  /v 1+b+c ; 


N+ij  1  , T  1  it  i  £  L  L— 1  ,  ,  L— £ 

V  /N+L-l N , 1+b+c  N  ,  v  / N+L-l N ,  be  \ 

l  ^  £  n  be  '  l  (  9 


£  1+b+c 


.  L  ....  N+L-l 

<I+sfe>  (1+^ 


.  t(l+b)(l+c)lW'L  1 
(bc)N_1(l+b+c)L 
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Therefore  finally, 


P  =  (i±k) 

D  V  b  ' 


N-l 


L  - 

,  1+b  N  1+b+c 

(T+b+^  6 


N-l 


Cl+c) 


N+L-l 


N-l 

I  < 

£«  0 


N+L-l 

L+£ 


)  e 


£+l'‘l+c; 


cd 

1+b+c 

e 


(^)£ 


e  /  b  cd  ' \ 
£+l4+c  1+b+c ;  I 


(5-17) 


The  corresponding  false  alarm  probability  is  still  given  by  (5-11),  and  it  can 
be  shown  directly  that  (5-17)  converges  to  (5-11)  as  b+0,  and  to  (4-22)  as 
d+0. 

We  forego  a  similar  analysis  of  the  Swerliug  case  3  problem  and  return  to 
(5-12)  with  the  assumption  that  M>N.  Then  (4-4)  can  be  used  for  all  the 
terras,  and  we  obtain 


(1+c) 


N+L-l 


N+L-l 

l 

£=0 


, N+L-l 
1  ; 


(  Hc- 

±+b+c 


M-N+£  U-N+£ 

)  l 

m=0 


ra 


m 


cd 


N+ra-£  1+b+c 


-) 


(5-18) 
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It  is  useful  to  rearrange  the  sums  in  (5-18)  in  order  to  get  an  increasing 
sequence  of  S-functions  (to  permit  their  recursive  computation),  hence  we 
write  N+m-£  =  n+1,  and  eliminate  m.  The  lower  limit  on  the  new  index,  n,  is 
set  by  the  S-functions  themselves  at  n=0,  and  the  upper  limit  is  seen  to  be 
M-l.  In  reordering  these  sums,  we  note  that  the  original  sums  over  m  and  l 
are  automatically  limited  by  their  binomial  coefficients,  since  (”)  is  zero  if 
either  k  or  n-k  is  negative.  Thus,  can  be  written 


PD  = 


M-l 

l  I 

n=0  <  £ 


,N+L-1. ,  M-N+£  .  l\  .  b  ,n 
(  £  )(n+l+£-N}  X  f 


( 


cd 


n+1  v 1+b+c 


) 


(5-19) 


where 


F  = 
o 


(1+c) 


M-N-L 


b^d+b+c)™ 


and 

-  he 

x  =  1+b+c  ’ 


The  lower  limit  of  the  £-sum  in  (5-19)  depends  upon  n. 
Next,  we  use  the  identity 


C-B 

x 


(1+x) 


A-C 


A+£ 

A+B-C 


) 


(5-20) 
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Ml  II  II  'I  iW  <l<lk>»l»K!t|llt>li||,H|llill  till!  'till  illH  lb  I  Hi  I|  I'M  u"l  1  II  ml  iN  I' ll  ll  III  I  III  1  'lb  •bliMbm  UlllMlI'lilll '.  | 


p  -  c^CL+b)^  MyN  (M-N  be  * 

D  "  a+b^)^-1  »-o  *  1+b*= 


C  ,N+L-mwl+b  “  .  cd  v 

x  i„  (  L+n  )(~ >  bn+l  T+b+c  ’ 
n=*0 


(5-21) 


In  Swerling  case  2  (M-N) ,  the  £-sum  In  (5-21)  reduces  a  single  term;  and 


,  c  (1+b) 

D  (l+b+c)N+L~1 


N  1  ,N+L-lwl+b  n 


V  cI'mj'J-v±I£)  s  (  Cd  ) 

L  '  T  +r.  «  '  VruK+b'  * 

n=0 


(5-22) 


The  false  alarm  probability  for  any  M  is  still  given  by  (5-11),  to  which 
(5-21)  reduces  when  b=0. 

For  Swerling  case  4,  (5-21)  must  be  used  with  M=2N.  It  should  be  noted 
that  (5-21)  can  be  obtained  somewhat  more  directly  by  averaging  (4-4)  over  the 
pdf, 

fo(y)  •  I  £l<c  -  d>  -  <5-23> 


for  the  threshold  y.  The  required  . .Legrals, 


/  Sf+1 (Xx)fL(x,d)  dx, 


are  evaluated  by  writing 
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fL(x,d)  -  dx  PL(x,d), 


then  integrating  by  parts,  and  using  (4-16). 

So  far,  we  nave  been  discussing  extensions  of  (5-8)  in  which  the  signal 
parameter,  a,  has  been  randomized.  In  the  dual  problem,  the  parameter  d  is 
randomized  instead.  This  represents  linear  CFAR  with  a  fixed  signal 
parameter,  a,  but  with  fluctuating  signal  components  in  the  randomized 
threshold.  For  d  we  assume  the  pdf 


fo(d) 


(5-24) 


and  average  (5-8)  with  respect  to  it.  The  result  is 


(1+c) 


N+L-l 


N+L-l  hj,  T  . 

I  (“t1) 

2-0  1 


lr  D  ,  a  .  ,1+c 

c  '  PN-t*y»l+c*  gK(" 


ch 


y) 


1+c 

ch 


dy  (5-25) 


after  the  change  of  variable  d+(l+c)y/c.  We  can  use  (4-16)  to  evaluate  the 
integral,  which  corresponds  to  ordinary  linear  CFAR  and  fixed  signal 
parameter,  with  the  variable  changes: 

N+N-A 

L+K 
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a  + 


a 

1+c 


c  + 


ch 

1+c 


Therefore 


(  r\  /  a  >  t  1+C  v  1+C  j 
J  PN-£(y’  1+?'  gK(-dh  y)  “ch  dy 


N+K-£-l  N+K-£-l  .  .  ,  m 

_  /  1+c  s  r  ,N+K-£-lwch 

-*■  M  1 _ i _ l_  ) 


'  1+c+ch ' 


l  <  v  >0  S, 

m=N-£  m  1+C  1 


( — - ) 


m  /vl+c/  m+-£+l-N  1+c+ch 


We  assume  that  K>L  to  assure  the  validity  of  this  formula  for  all  £(in 
particular,  £=N+L-1).  When  (4-16)  was  derived,  the  parameter  there  called  N 
was  non-negative,  and  we  are  now  applying  to  a  more  general  case.  A  review  of 
the  derivation  of  (4-16)  will  show  that  it  is  valid  as  long  as  N+L-1>0,  which 
explains  the  requirement  K>L.  This  excludes  Swerling  case  1  type  models  for 
the  fluctuating  signal  component  in  the  threshold. 

Substituting  in  (5-25)  we  get 


/ -i  ,  _  L  N+L— 1  -T  T  ,  .  ,  .  ,  £ 

r.  _  i  (1+c)  r  ,N+L-lNr  1+c+ch, 

p° ~ 1 ' Jo  (  *  )tc~1  * 


1  ,N+K-£-lu  ch  .m  ,  a  . 

X  (  m  ^  1+ch  Sm+£+l-N  ^l+c-;-'^ 
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The  sum  on  m  has  been  allowed  to  start  at  m=0,  since  the  S-function  will 
vanish  in  the  extra  terms.  With  substitution  £-»N+L-l-£,  this  becomes 


PD  "  1  - 


N+L-l  K-L  N+L-l  , 

,  1+c  v  Y  / N+L-l v  ,1  1+c  , 

^  i  £  >  Lc  1+c+chJ 


'  1+c+ch ^ 


K-L+£  .  ,  m 

V  ,K-L+4v  ,ch»  g  / _ 

Ln  K  m  ' ' 1+c '  L+m-£  ' 1+c+ch 

m=0 


) 


(5-26) 


Equation  (5-26)  can  be  obtained  from  (5-] 8)  by  reversing  the  roles  of  the 
original  random  variables,  complimenting  the  probability,  and  making  the 
parameter  substitutions 


N+L 

L+N 

ll+K 


b-*h 

d+a 

c+l/c 


(5-27) 


A  form  analogous  to  ’5-21)  can  be  obtained  from  (5-26),  either  by  a  parallel 
calculation  or  by  making  the  variable  substitutions  (5-29)  in  (5-21).  Either 
way,  the  result  is 


1  - 


K/i.u,\N+L-K 

c  (1+h) 


(1+c+ch) 


N+L-l 


K— L  ,r  .  .  SL 

i  c';lki  h 


£=0 


1+c+ch 


) 


x 


L-l+.l 


n=0 


.N+L-l+JL  ,  ,.n  c 

(  s+n  X'*10  s 


< 


n+1  1+c+ch 


) 


(5-28) 
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It  can  be  seen  that  (5-28)  reduces  to  (4-16)  when  h=0,  since  then  only  the 
tern  £=0  survives. 


The  false  alarm  probability  is  obtained  from  (5-28)  by  putting  a=0,  which 
has  only  the  effect  of  replacing  the  S-function  by  unity.  But  then. 


L-l+£  NfL-l+£ 

i  <^;m)<c-*h)” ,  <c4ch)-»  i  (‘,+L:i+i)(C^h)" 

n=0  U"  n-N 


(1+c+ch) 

(c+ch)f 


N+L-l+£ 


N-l 

1  r  ,N+L-1+£W  ,  ,.n 

- s  l  (  )(c+ch)  , 

(c+ch)  n=0 


(5-29) 


I  (K:L)h"=  (l+h)K'L  . 
£='0  * 


The  result  is  that  the  first  term  of  (5-29)  eventually  cancels  the  unity  in 
(5-28)  and  the  false  alarm  probability  simplifies  to 


/,  ,u\Ij— K  K— L  £ 

p  sx  — LL ill -  y  (K  — )  x 

FA  N+T.-1  4  v  5  '  a 


FA  ,  SN+L-1  on  i  1+c+ch 

(1+c+ch)  £=»0 


r  ,N+L-l+£v,  ,^n 

x  )  (  )(c+ch) 


(5-30) 
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This  expression  reduces,  of  course,  to  (4-17)  when  h=*0. 

The  final  situation,  in  which  both  of  the  signal  parameters,  a  and  d, 
are  randomized  is  analyzed  first  by  averaging  (5-28)  with  respect  to  the  usual 
pdf  for  a: 


£o(a>  =  k  %<!> 


In  this  way,  a  formula  valid  for  all  M  is  obtained,  although  the  constraint 
K>L  is  still  in  force.  Only  the  S-function  in  (5-28)  depends  on  a,  and  the 
use  of  (4-20)  gives  us 


/  Sn+l(l+c+ch)  %^b^  d(b) 


,  1+c-tch  >M  r  b 

^1+b+c+ctv  m  1+b+c+clr 

m=0 


(5-31) 


When  (5-31)  is  used  to  average  (5-28),  the  result  is 


nN+l~K  ,, ,  ,  ,  -M-N-L+l 
D  _  ,  c  (1+h)  (1+c+ch) 

D  M 

(1+b+c+ch) 


K~L  v_  T  U  l  L-l+£ 

l  J0  - 


69 


X 


m=0 


,11-1+nK  >  b  ' 
m  '^1+b+c+ch' 


m 


(5-32) 


and  the  probability  of  false  alarm  is  still  given  by  (5-30). 

To  reverse  the  constraints,  we  can  obtain  a  formula  valid  for  all  K  and 
11  >N  by  averaging  (5-21)  with  respect  to  the  pdf  of  d,  namely  (5-24).  The 
resulting  detection  formula  is 
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’D 


N-l 
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(l+b)N^Wc)K~N'L+1 

(1+b+c+ch)^ 
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x 


11-N 


l 

£=0 


)( 


be 


1+b+c 
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N-l+£ 
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n=0 


X 


X 


n 


I 

m=0 


•K-l+nu ,  ch 
m  ^  1+b+c+ch 


m 

) 


(5-33) 


Expression  (5-33)  reduces  to  (4-26)  when  h=0,  as  does  (5-21)  in  the 
corresponding  limit,  d=0.  It  is  possible  to  show  that  (5-32)  and  (5-33)  are 
equivalent  when  both  K>L  and  11>N. 

The  false  alarm  probability  corresponds  to  b=0: 
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X 


cN-1(i+c)K-N-L+1  N-l  n 

--  -K  L  K  L+n  ncJ 


(1+c+ch) 


fK-l+mv  ch  u 
k  m  nl+c+ch; 


(5-34) 


We  did  not  work  out  the  detection  and  false  alarm  probabilities  for  the 
analogous  fixed-signal  case,  with  K<L,  but  (5-34)  is  identical  to  the  result 
of  averaging  (5-11)  over  the  pdf  (5-24)  for  d. 

Formulas  (5-32)  and  (5-33)  generalize  some  results  of  Rickard  and 
Dillard/11^  and  also  Wishner’s  analysis of  the  normalized  periodogram 
detector. 

VI*  RELATION  TO  OTHER  FORMS 

A  striking  feature  of  the  finite  sum  detection  formulas  is  the  diversity 
of  equivalent  forms  they  can  assume.  Solutions  to  the  same  problem,  such  as 
(4-21)  and  (4-26),  deduced  from  different  starting  points  can  be  reconciled 
only  with  difficulty.  In  this  section,  we  relate  a  few  of  our  results  to 
quite  different  formulas  which  have  appeared  in  the  literature. 

Consider  Eq.  (4-16)  which  describes  linear  CFAR  with  a  fixed  signal 
parameter,  a.  It  is 


1  -  P  = - - - 

D  a+c)S+L-) 


which  can  be  written 


,N+L-1,  £ 

N+£  )C  S£+l  ^1+c5  » 


(6-1) 
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N 

•) 


(1 


V 


FNa-D 


L-l 

l  T  (L-l,  £) 
£=0 


(6-2) 


where 


Tn(L,£) 


1  ,N+L.  2  C 

,L  C  SJ1+I1+C) 

(l+c; 


(6-3) 


How  fn(L)  can  be  rewritten  as  a  sum  of  terms,  each  of  which  depends  only  on  £, 
and  not  also  on  L,  as  do  the  terms  T^(L,  £).  To  accomplish  this  we  write 

F»‘L>  ■  X  tn(l>«  ■  VL-U  +  VL-« 

£=0  £=0 


L-l 

-  T  (L,L)  +  l  [T  (L,£)  -  T  (L-l,£)]  +  F„(L-1), 
N  £=0  N  N  N 


or 


Fn(L)  =  AN(L)  +  FN(L-1), 


(6-4) 


where 


L-l 


AfL)  =  T„(L,L)  +  l  [T„(L,£)  -  T„(L-1,£)3 
h  N  £=0  ” 
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(6-5) 


ill 


for  L>1.  Iterating  (6-4)  we  obtain 


tN(L)  =  AN(L)  +  an(L_1}+***+AN(1>  +  V0)  * 


Defining 


V0)  =  Fn(0)  =  IKt(0,0)  , 


(6-6) 


we  can  write 


VL> 


(6-7) 


which  has  the  desired  characteristic. 
As  applied  to  (6—2),  we  have 


1 

(l+e)L 


,N+L.  *  c  ,  N 


(6-8) 


where 


One  finds  that 
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-u  r  ,N4L-1.  (cu) * 
i=o 


(6-9) 


which  holds  alsc  for  L=0.  Tnerefore 


F„(L)  =  e  U  l 


r  (cu)m 

i  Vfi+T- -1  1  ml 


£=0  (14c) 1  ra=0  llHr'1  m! 


and  L  now  appears  only  as  the  upper  limit  of  the  outer  sum. 


me  finite  nt-sum  is  a  special  case  of  the  confluent  hypergeonetric 

/ 12\ 

function;  in  fact  it  is  a  Laguerre  polynomial.  The  basic  definition  is 


L(k)(x)  =  *  (k4n} 

m  lx)  i  'k+a^  n! 


(6-10) 


hence 


c  n  -v  _  ~u  r  1  -  (N-l) ,  , 

fn(L)  -  e  l  - -£  L£  (-cu)  , 


£=0  (14c) 


and  finally 


„  K  "iXT  L~l 

**-!-&>  *  ^  l 


£=0  (14c) 


,(N-1)  ,  ca  . 
t  14c } 


(6-11) 


The  Laguerre  polynomials  can  be  computed  from  the  second-order  recursion 

.  (12) 

l elation 


L<k)(-x)  =  (2  +  L(£k>(-x)  ~  (1  +  L(£k>  (-x) 


(6-12) 


together  with 

L^k)(-x)  =  1 

(6-13) 

Ljk)(-x)  »  k+l+x  . 


(13) 

The  particular  case,  N=l,  of  formula  (6-11)  was  given  by  Finn,  whose 
L£(x)  equals  it  L£^(-x),  in  our  notation. 

It  is  interesting  that  the  probability  of  false  alarm  obtained  from 
(6-11)  appears  in  the  fora 


P  =  1  - 
FA  vl+c' 


L-l 

r 

L 

£=  0 


<T> 


l 

£  * 

(1+c) 


(6-14) 


which  can  be  transformed  into  (4-17)  by  the  use  of  (4-27). 

Formula  (4-4)  can  be  treated  the  same  way,  when  M>'J.  According  to  (4-4) 


Pn  =  K  cm) 


hence 
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*  ,„x  uu  -u  (K-l) !  ,  (N)  ,  N 

VK)  (l+b)K  6  (N+K-l ) !  LK-1  ("bu) 


Finally, 
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=  SN(u)  +  bu  e  Lk-1  ('bu)  > 


M_y_1  k!  1  (N)  by_ 

k=0  (N+k)!  (l+b)k  k  l+b 

(6-16) 

This  time,  the  false  alarm  probability  (b=0)  is  simply 

PFA  '  Vy)  > 

the  same  as  (4-5). 

We  demonstrate  another  kind  of  equivalence  now,  reducing  an  infinite 
series  solution  for  the  fixed-signal,  linear  CFAR  problem  to  our  finite  sum 


and  therefore 


„  s  (-Z-)  +  —  (-Z-) 

Vl+b;  l+b  'l+b'' 


_  _I_ 

l+b 
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(4-16).  The  former  solution  is  obtained  by  averaging  Fehlners'  series: 


PN(y»a) 


(6-17) 


over  the  pdf 


£„<y> 


fcr  the  threshold.  Series  (6-17)  can  be  obtained  from  (3-7)  by  expanding  the 
factor  exp(a/t)  in  a  power  series  and  integrating  term  by  term.  These 
integrals  are  performed  by  using  (3-14). 

According  to  (4-20), 
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(6-18) 
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The  last  step  in  the  derivation  of  (6-18)  is  an  application  of  (4-27). 
Another  transformation  yields 
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(6-19) 


When  (6-19)  is  used  to  average  (6-17),  the  result  is 
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^  N+m  ^  iF/N+L‘*  N+m+l;  1+c)  » 


in  terms  of  the  confluent  hypergeometric  function/15^  Applying  Kummer's 
first  transformation/16^ 
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The  new  hypergeometric  function  is  simply  a  polynomial,  since  its  first 
argument,  m+l-L,  is  a  non-positive  integer.  Altogether, 
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(6-21) 


Next,  we  replace  is  by  the  new  index  n  =  m+4  ,  and  find 
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which  is  -i&).  Similarly  by  averaging  (6-17)  ever  a,  u-ing  the  usual  pdf 
with  M>N,  (4-4)  can  be  derived. 
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APPENDIX  1 


The  binomial  identity 


n 


l 

m=0 


,M-l+nK  m 

(  „  >  x 


(l-x)n 


(A-l) 


is  proved  by  taking  the  factor  (l-x)n  inside  the  sum  and  expanding  again,  as 
follows : 
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£=0  f  m=0 


(M+n)(?j")(-i)£-In 

m  £~m 


In  the  last  step,  we  have  rearranged  the  sum  in  increasing  powers  of  x.  now 
the  m-sum  is 
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M-l+£, 

■  0  j 


(A-2) 


according  to  a  standard  binomial  identity,  which  completes  the  proof.  The 
standard  identity  used  here  may  be  written  in  the  form 


l  (g)(£"^_S)(-l)S  =  (-lAk”*>  ,  (A-3) 

s=0 


which  is  valid  when  k>£,  and  is  proved  by  use  of  the  obvious  identity 


(l+x)k  (l+x)“£  =  (l+x)k  1  . 

Each  factor  on  the  left  is  expanded  in  a  binomial  series  and  the  resul 
double  sum  reordered  in  powers  of  x.  The  desired  result  follows  by  equating 
to  the  expansion  of  the  right  side. 

To  prove  (5-20),  we  need  an  intermediate  result.  Consider  the  identity 


m 

[’+x(l+y}]m  =  (!■*)“  (1  +  yg-) 


(A-4) 


and  expand  the  left  side  and  the  right-most  factor  on  the  right  side: 
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I  o  x  (1+y >  -  u-*)m  l 


r  / m .  .  .  n  ,,  ,  m-n 

=  L  <.nKxy;  (1+x) 

n 


Explicit,  summation  limits  are  not  used  because  the  sums  are  limited 
automatically  by  the  binomial  coefficients,  as  mentioned  in  connection  with 

g 

Eq.  (5-19).  Now  we  multiply  both  sides  by  (1+y)  and  expand  some  more.  Th 
left  side  becomes 
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and  the  right  side  becomes 


/i .  \S  r  ,m..  n,,  Nm-n  n 

(1+y)  l  (n)  x  (1+x)  y 
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-  I  (“)  *"  7  (s)  yn+j 


we  introduce  the  new  index,  k=n+j,  and  eliminate  n,  so  that  the  last  double 
sum  expression  becomes 
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Equating  powers  of  y,  we  obtain  the  desired  result: 
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It  is  necessary  only  that  m,  s  and  k  be  non-negative,  and  the  ranges  of  the 
sums  over  £  aid  j,  as  limited  by  the  binomial  coefficients,  will  depend  on  the 
relative  magnitudes  of  m,  s  and  k. 

Now  consider  (5-20): 
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Using  (A-5)  of  che  left  side  we  get 


xCd+x)A"C 


I  («)(CA  )(^) 

3  J 


j 


and  applied  to  the  right: 
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In  the  last  step  we  replaced  j  by  A-C+j,  and  this  completes  the  proof  since 
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